Bosons with attractive interactions in a trap: Is the ground state fragmented ? 
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Possible fragmentation of a Bose-Einstein condensate with negative scattering length is investi- 
gated using a simple two-level model. Our results indicate that fragmentation does not take place 
for values of the coupling for which the system is metastable. We also comment on the possibility 
of realizing a fragmented condensate in trapping potentials other than an harmonic one. 
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Bose-Einstein condensation in dilute atomic gases has 
been an area of intense activity since the first experi- 
mental realizations of these systems |l] ||]. The experi- 
ments have so far been carried out with the alkali atoms 
87 Rb, 23 Na, and 7 Li. While the s-wave scattering length 
is positive in the first two cases, it is negative for Li, im- 
plying that the effective atom-atom interaction is attrac- 
tive. A homogeneous system of these atoms will collapse 
to a dense state before densities where a Bose-Einstein 
condensate can form are reached However, in the 
experiments Bose condensation is realized under inho- 
mogeneous conditions, and theoretical calculations show 
that the system may then be metastable as long as the 
number of condensed particles is below a critical value, 
N c ~ 10 3 for 7 Li for the trap used in the experi- 
ments. The latest experiments are consistent with this 
result 0. 

The homogeneous Bose gas with attractive interactions 
was considered in a paper by Nozieres and Saint James 
§ . While their main interest was the onset of a BCS-like 
transition, they also asked whether a fragmented Bose 
condensate would form. In the homogeneous case, frag- 
mentation means that the particle distribution in mo- 
mentum space, rik has a sharp peak near k = which 
extends over a number of (k, — k) states large compared 
to 1, but small compared to the particle number N: on 
a macroscopic scale, it looks like a 5-function. That this 
is a reasonable question can be seen within the Hartree- 
Fock approximation: for bosons with an attractive inter- 
action the Fock term makes it energetically favorable to 
spread the particles over several states. Still, the authors 
of Ref. H found that fragmentation of the condensate 
does not take place in a homogeneous system. However, 
the physics of the inhomogeneous Bose gas is different, 
and it seems worthwhile to ask the same question again 
in this context. 

We want to study the ground state of bosons with 
a negative s-wave scattering length in an isotropic har- 
monic oscillator potential. Since the question we are 
addressing is whether the condensate can fragment, we 
must allow for (at least) two single-particle states in the 
calculation. We will consider the two normalized single- 
particle wave functions 
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where r 2 = x 2 + y 2 + z 2 and b is a variational parame- 
ter. These have the same shape as the ground state and 
the first excited state with angular momentum projec- 
tion mi = of the harmonic oscillator well. There will 
be interactions between particles in the different levels 
and also between particles in the same level. For fixed b 
the many-body Hamiltonian is 
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where i,j,k,l = 0, 1, t is the kinetic energy operator, 

(4) 



where a is the scattering length (assumed negative) and 
m is the boson mass. We will in the following solve this 
two-level problem for arbitrary values of 6, thus obtain- 
ing the ground state energy of the model as a function 
of this variational parameter. In the next step b is deter- 
mined by minimizing the energy. This simple two-level 
system can be solved exactly and is in fact a special case 
of the so-called Lipkin model Q which has been widely 
used in nuclear physics as a tool for studying many-body 
approximation methods. 

We now proceed by writing out the Hamiltonian in Eq. 
(^) in more detail. The single-particle matrix elements 
are given by 
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where x = b/a osc , and fl^ sc = h/mu). Since ipo and ipi are 
real and of opposite parity, the only non-zero two-particle 
matrix elements are 
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By defining 2e = Tiuj{x + 1/x ), we obtain 
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which is of the same form as the Hamiltonian for the 
Lipkin model Introducing the so-called quasi-spin 
operators 
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which obey the SU(2) algebra of ordinary spin, that is 
[J+, J_] = 2J 2 and [J z , J±] = ± J±, the Hamiltonian can 
be written as 
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where N = a ao + a\ai is the particle number operator. 
The dimensionless parameter in the problem is seen to 
be N\Uo\/e. Furthermore, we introduce the "Cartesian" 
components of the quasi-spin through J± = J x ±iJ y and 
rewrite H as 



H = eJ z + 2eN--\U \ 
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Since the total spin j 2 = (j + j_ + j_ J+)/2 + j 2 com- 
mutes with iJ, the eigenstates and eigenvalues can be 
obtained exactly by diagonalizing (2 J+ 1) x (2J + 1) ma- 
trices. For N particles, the ground state is found among 
configurations having J = N/2. We can also solve the 
problem analytically in the semi-classical approximation. 
In this approach the angular momentum operators be- 
come c-numbers and we take 
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Upon substituting these expressions in H, and taking 
N/2(N/2 - 1) w 7V 2 /4, N - 1 w AT, we find the semi- 
classical expression for the energy of the system 
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This quantity is for fixed 9 clearly minimized by 4> = 0, 
and thus the energy per particle can be written as 



E 1 
N = 2 [ ' 



N\U \ 



5N\U \ 



32 



-T) 



1 /19.2V 



2 V 16 



1 \U Q \ + 2e 



(21) 



with r\ = cos 9. Equilibrium occurs at the value 770 of 77 
given by 
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and the equilibrium energy is found to be 

_i[-3e + (j\r-i)|E7b|], mm <2 , 

^hl^+j^ + (QN-5)\U \l ^>2. 
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For the population of the lowest single-particle state ifjo 
we get 
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Thus we see that the extent to which the condensate is 
fragmented is determined by the ratio 
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The quantity N\a\/a osc will from here on be called the 
effective coupling. A large effective coupling will give a 
correspondingly large population of the state ip±. How- 
ever, since the two-body interaction is attractive the sys- 
tem will collapse if the effective coupling exceeds a crit- 
ical value. From a variational calculation with only the 
state tpo taken into account, the critical value is found 
to be N\a\/a osc « 0.67 in reasonable agreement 

with the stability criterion N]a\/a osc < 0.58 obtained by 
Ruprecht et al. from numerical solutions of the time- 
dependent Gross-Pitaevskii equation 0. In Fig. [I] we 
show results for both the exact and the semi-classical cal- 
culation of the ground state energy of the Hamiltonian 
in Eq. Jl5|). The exact solution was obtained numeri- 
cally. The quantity H/N depends on TV essentially only 
through the effective coupling N\a\/a osc , and we checked 
that it was sufficient to do the calculations with TV = 50 
particles. This made numerical calculations easy, as we 
simply had to diagonalize 51 x 51 matrices for various val- 
ues of the variational parameter x — b/a osc . From Fig. 
|l| it is seen that the semi-classical solution is very close 
to the exact one. The critical value of the effective cou- 
pling was again found to be rj 0.67. At this low effective 
coupling, the interaction is not strong enough to excite 
particles into the state ipi to any appreciable extent, and 
the result is therefore nearly identical to the calculation 
of Ref . . This is reflected in the occupation number 
for ipo, shown as a function of N\a\/a osc in Fig. ^. For 
all effective couplings which allow for a metastable Bose 
condensate, all particles condense in i/jq. The depletion 
reaches ~ 10% at N\a\/a osc m 3.5, and approaches 40% 
in the limit N\a\/a osc 3> 1. The derivative of the semi- 
classical result for TV /TV is discontinuous at JV|J7o|/e = 2, 
while from Fig. ||this is not the case for the exact result. 
This is due to the finite value of TV in the exact calcula- 
tion, as we checked that the edge at N\Uo\/e = 2 became 
sharper when we increased TV. 

Our choice of single-particle wave functions led to spe- 
cific ratios of the interaction strengths Voooo> Vim, and 
Viioo- As one could imagine alternative choices for the 
wave functions (experimentally one could alter the trap- 
ping potential), it is worthwhile to check how a change 
in these ratios affects our results. We use the same form 
of the Hamiltonian in Eq. (11), but take 
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and 
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In the semi-classical approximation the energy per par- 
ticle becomes 
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neglecting constant terms of order 1 compared with TV. 
Minimizing with respect to r\ we find that the population 
of the lowest single-particle state is given by 
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and we see that we could get a significant fragmentation 
of the condensate even with N\Uo\/e ~ 1 if we could 
make (3 large. However, the elementary inequality 
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which gives 1 + a > 2/3, will in practical situations pro- 
hibit us from designing the two-body matrix elements so 
that significant fragmentation is obtained. If 1 + a = 2/3 
then \ipo\ 2 — IVM 2 ! i- e - "00 an d V'l are identical and there 
is no fragmentation. 

To conclude, we have examined a simple two-level 
model to see whether an attractive interaction favors 
the formation of a fragmented Bose-Einstcin condensate. 
Our results indicate that the answer to this question is 
no: at the coupling strengths where a metastable state 
exists, all particles condense in the state rpQ. We have 
only considered the case of condensation in an isotropic 
harmonic oscillator potential, and the situation could be 
different if one were free to design the relative strengths 
of the two-body matrix elements involved in the calcu- 
lation. However, the inequality ([H]) must be satisfied, 
and this seems to prevent us from choosing the strengths 
so that fragmentation is obtained. Thus the conclusion 
remains the same as in the homogeneous case considered 
by Nozieres and St. James ||: the system will collapse 
before a fragmented Bose condensate can form. We re- 
mark that our results are consistent with those of Wilkin 
et al. p^ |, who found no fragmentation for bosons with 
weakly attractive interactions in the absence of rotation. 

0E is grateful for the hospitality of Nordita, where this 
work was carried out. 
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FIG. 1. Comparison of the exact and semi-classical result for the energy per particle. 
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FIG. 2. Population of the lowest single-particle level as a function of the effective coupling. The exact result is for iV = 50 
particles. 



•5 



